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ABSTRACT 


Au/n-Si Schottky barrier diodes have been fabricated by depositing a gold film 
of Ri 600 A thickness, using a thermal evaporation technique in vacuum « 10“® torr, on to 
a cleaned n-type silicon wafer of (111) orientation Their current- voltage (TV) character- 
istics have been measured at room temperature and analysed on the basis of thermionic 
emission-diffusion (TED) mechanism, using a computer program The parameters, namely 
barrier height, ideality factor and series resistance are shown to vary from diode to diode, 
the reasonable values being 4>b ~ 0.79 V, 17 Ri 1.17 and Rs « 200 — 300 fi. 

The values of barrier height obtained from the Capacitance- Voltage (C-V) character- 
istics are invariably higher than those derived from the (TV) data. This discrepancy and 
high values of ideality factor are attributed to the prevailing barrier inhomogeneities at 
the gold-silicon interface. Based on X-ray diffraction, scanning electron micrographs and 
Rutherford back scattering spectrum, it is shown that gold does not react with silicon 
and thus form a pure metal-semiconductor junction. Further annealing of the diodes at 
450° C for 30 minutes leads to the formation of gold globules and deterioration in their 
characteristics. 

TV characteristics of a Schottky diode have also been simulated using a TED cur- 
rent expression and assuming a Gaussian distribution of barrier heights. Finally^ it is 
shown that the decrease of barrier height, abnormal increase of ideality factor and non- 
linearity in the activation energy plot of the saturation current (usually observed in prac- 



tice) can be explained by assuming the presence of barrier height inhomogeneities at the 
metal-semiconductor interface. 
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Chapter 1 


Introduction 


A rectifying metal-semiconductor (M-S) contact is known as a Schottky barrier diode and it 
exhibits asymmetrical electrical characteristics [1-2]. During the last century these contacts 
have experienced three major periods of industrial importance; first as radio detectors in 
early 1900s, then as radar detectors during the second world war and later on as clamping 
diodes, microwave diodes and gates of microwave transistors [3-5]. The first M-S device 
was the point contact diode made by pressing a metal whisker to a semiconductor crystal 
piece. These days a high vacuum system is employed to fabricate reproducible planar 
metal-semiconductor interface with uniform contact potential and current distribution over 
the junction area. 


1.1 Formation of Schottky barrier 

When a metal comes in intimate contact with a semiconductor, a potential barrier is created 
by separation of charges at the M-S interface. Concurrently, a high resistance region devoid 
of mobile charges is formed within the semiconductor near the interface Schottky and 
Mott explained the mechanism of barrier formation on the basis of the difference in the 
work functions of the metal and the semiconductor. Fig. 1.1a illustrates the energy band 
diagram of a case in which the metal has a larger work function compared to (say) the n- 
type semiconductor. When the metal and n-type semiconductor are brought into contact. 
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electrons transfer take place from semiconductor to metal until Fermi levels get aligned 
and equilibrium is established. Electrons movement towards metal causes depletion in the 
semiconductor near the interface. Since the separation between the conduction band edge 
Ec and Fermi level Ep increases with decreasing electron concentration, the conduction 
band edge turns upwards near the interface. The migration of electrons leaves positively 
charged donors in the semiconductor and forms a thin sheet of negative charges at the metal 
surface. This produces an electric field in the direction opposite to electrons flow, i.e., from 
the semiconductor to the metal The width of the space charge layer in the semiconductor 
is appreciable because donor density being several orders of magnitude smaller than the 
electron concentration in the metal. Since the band gap of the semiconductor remains 
unchanged after the contact, the valence band edge Ey also moves upwards identically to 
the conduction band edge Ec- Further, the vacuum level m the semiconductor follows a 
similar trend and bends exactly as Ec- Thus, for a M-S system in thermal equihbrilium 
the vacuum level remains continuous across the transition region. The amount of band 
bending is then just equal to the difference between the two vacuum levels ( see, eg., 
Fig.l lb). Thus, 

4^s 

where K is the built in potential at the junction, (pm and ps are the work functions (in 
units of eV ) of metal and semiconductor, respectively. The Vi represents the potential 
barrier which an electron has to overcome while moving from the semiconductor to the 
metal. The corresponding energy is qVi. The barrier (pj,) as seen from the metal side is 
given by p^ = p^- Xs ox = pm- {Ps - Pn), since p, = Xs + </>n»Using eq. (1.1), we have 

Pb = qV^ + Pn (1 2 ) 

where Xs is the electron affinity and pn — Ec — Ep represents the penetration of the Fermi 
level in the band gap of the semiconductor. Obviously, pb > qV.^. The exact shape of 
the barrier is determined from the charge (or donor ) distribution within the space charge 
layer. Schottky assumed the semiconductor to be uniformly doped, giving a constant 
charge density in the depletion region. The electric field strength therefore rises linearly 
with distance from the edge of the space charge layer. Mott assumed a thin layer devoid of 
charges sandwiched between a uniformly doped semiconductor and the metal. The electric 
field strength within the so called depletion region is constant. Fig. 1.1 la shows the thermal 



METAL 



(0 ) 
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Fig. 1.1 Electron energy band diagrams of metal and n-type semiconductor having <pm > <Ps 
(a) neutral metal and semiconductor separated from each other, (b) after the contact and 
attaining thermal equilibrium . 
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equilibrium energy band diagram of the M-S contact The depletion region of width w has 
few mobile carri^ and therefore displays very high resistance. Consequently a large part 
of the externally applied voltage drops across the depletion region. If a negative voltage 
(V/r) is applied at the semiconductor, the depletion region width decreases and voltage 
across it falls from K to (K - Vp)- The electrons now see a reduced barrier and move 
from the semiconductor towards the metal causing a current flow. Since no voltage drop 
occurs within the metal, remains unaffected and there is no electron flow from metal 
to semiconductor. Thus, M-S junction exhibits rectifying properties 


1.2 The concept of surface states 

As per eq.(l.l), the barrier height should increase linearly with but, is found to be not 

so in practice. The discrepancy has been explained by Bardeen on the basis of suface 
states. Accordingly, broken covalent (or dangling) bonds at the semiconductor surface give 
rise to localized energy states, called surface states. These are believed to be distributed in 
the band gap and characterized by a neutral level (po- That means when these is no band 
bending in the semiconductor, the surface states are occupied right up to (po making the 
surface neutral, i.e., a flat band condition shown in Fig. 1.2a. The states below (po are donor 
like and those above po behave like acceptors. Their presence leads to the redistribution of 
the charge near the semiconductor surface and affect the barrier height. The equilibrium 
is established when sufficient number of electrons from the semiconductor adjacent to the 
surface occupy the localized states lying above po- The surface then becomes negatively 
charged and a depletion layer consisting of ionized donors is created near the semiconductor 
surface. Thus, a potential is formed even in the absence of the metal contact as depicted 
in Fig 1.2b. If a metal is now brought into contact with the semiconductor and equilibrium 
reached, the Fermi level in the semiconductor changes by an amount equal to the contact 
potential difference through transfer of electrons to the metal. If the density of surface 
states is high, the charge exchange takes place largely between the metal and the surface 
states and the space charge in the semiconductor nearly remains unchanged. As a result, 
barrier height becomes independent of the metal work function and is given by [6] 

P, = (Eg - Po) (1.3) 
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Fig. 1.2 Electron energy band diagrams of n-type semiconductor with surface (a) un- 
der flat band , (b) m thermal equilibrium with the bulk , and (c) in contact with a metal. 



In this case, the barrier height is said to be ’’pinned” by surface states. Eq. (1.3) is known 
as the Bardeen limit. 

1.3 Conduction mechanisms 

The current flows across a Schottky barrier diode because of charge transport from the 
semiconductor to the metal or vice-versa through several different mechanisms namely, 
thermionic emission-diff'usion, tunneling, generation or recombination, and minoriy carrier 
injection. In addition, there may be leakage current due to high electric held at the contact 
periphery or interface current due to traps at the metal-semiconductor interface. 

1.3.1 Thermionic emission-diffusion mechanism 

It is the dominant transport mechanism in the Schottky barrier diodes and involves drift 
and diffusion both [7]. The electrons in the conduction band are assumed to have kinetic 
energy only and barrier height (j)b is much larger than (kT/q). The net current due to flow 
of electrons from semiconductor to metal by overcoming the potential barrier of height 
at a forward bias V is given by [8,9] 



Where A** is the effective Richardson coi^ant, Ad is the diode area, T is the temperature 
in kelvin, k is the Boltzmann constant, q is the electronic charge and h is the Plancl^s 
constant For free electrons, the Richardson constant A** is 1.12 xlfl4 In deriving 

eq.(1.4) diffusion of carries is also included. However, the barrier height invariably increases 
with the forward bias. If the barrier height is assumed to vary linearly with bias, we have 

UV) = <l>bo + lV (1.5) 

Where 4>bo is the barrier height at zero bias and 7(= d4>i,/dV) is positive. Substituting the 
value of (t)b{V) into eq.(1.4) yields ' 



where 


Is^A^A**T^exp 


^ —<l<PbO \ 
^ kT J 


(1.7) 


and 1/77 = 1 — 7 . The parameter rj is called the ideality factor and usually has a value 
greater than unity. The neutral region of the semiconductor offers resistance (Rs) and so 
significant voltage drop occurs there at large forward currents. As a consequence, voltage 
across the Schottky barrier region is reduced by a factor IRg and resulting current becomes 


I = Is exp 


' q{V-IRs) 

' rjkT 


exp (- 


Now a plot of In (I) vs V deviates from a straight line and exhibits saturation at high 
forward voltages. 


1.3.2 Tunneling through the barrier 


This mechanism prevails when the semiconductor is heavily doped such that the Fermi level 
lies above the bottom of the conduction band and the depletion region width is very small. 
Consequently, at low temperatures, electrons having energy nearly equal to the Fermi level 
can tunnel from the semiconductor into the metal. The tunneling current through the 
barrier is of the form [9] 

It = Jtoexp(^)[l - ea:p(-^)] (1.9) 

where 



where m*(= m^mo) is the effective mass of electron and €s(= e^eo) is the permittivity of 
the semiconductor, nio is the electron rest mass and Nd is the donor concentration in m”®. 
The process is called field emission (FE) when Eoo » kT. At high temperatures where 
Eoo kT, the tunneling of thermally excited electrons occurs and the process is called 


thermionic-field emission (TFE). It may be noted that TFE also causes a reduction in the 
barrier height by [9] 

f3\^ 2 i 

= i-\ EioVi (1.12) 

Where Vd stands for the voltage corresponding to the band bending. Eq.(1.9) suggests that 
the lnlt/[l — exp(^^^)] versus V plot should yield a straight line with slope giving Eq and 
the intercept at zero bias the saturation current ho- 


1.3.3 Carrier generation and recombination 


At zero bias, the rate of electron-hole pair generation in the depletion region is balanced by 
the rate of recombination and therefore there exists a thermal equilibrium. The electron- 
hole product (np) is then equal to n^. However, with the application of the bias, np 
product departs from and a net generation or recombination of carriers occur depending 
upon the polarity. In case of a Schottky diode made on n-type silicon, substrate is 
forward biased the electrons are injected into the depletion region from the neutral bulk 
semiconductor and holes swept away from the metal. These additional electron-hole pairs 
combine in the depletion region to give a forward combination current. Such a process 
occurs most effectively through recombination centres (or traps) with energy levels lying 
near the centre of the band gap. The recombination current can then be described by [9] 

I, = U^xp(A) _ 1] (1,13) 

with 


_ qn^AdW 
2t 


(1.14) 


and 


n, = ^NcNvexp{-^) (1.15) 

Here w is the width of the depletion region, r is the carrier effective life time within the 
depletion region, rij is the intrinsic carrier concentration. Eg is the energy band gap, and 
Nc and Ny are the effective conduction and valence-band density of states, respectively. 
When the Schottky diode is reverse biased additional electron-hole pairs are generated in 
the depletion region and cause a generation current in the reverse direction. 


1.3.4 Minority carrier injection 


The Schottky barrier diode is a majority-carrier device under the low bias conditions. 
However, at large forward bias, minority carrier current contribution increases and become 
more than the diffusion current. For example when the barrier height of the semiconductor 
(say n-type) is more than half of its band gap the interface region adjacent to the metal 
becomes essentially p-type. Under a forward bias, the electrons from the semiconductor 
flow into the metal and some of the holes are swept into the neutral region. Thus, injection 
of holes occur from the metal into the semiconductor. It may be noted that the process 
IS equivalent to flow of electrons from the semiconductor valence band into the metal 
directly. The current Ip caused by hole injection from metal into the neutral region of 
n-type semiconductor of donor concentration Nd is given by [10] 


Ip — 


qAdPpnl 

N^Lp 


exp(^) 


(1.16) 


Where Dp is the hole diffusion constant and Lp = \jDpTp is the hole diffusion length, 
represents the hole lifetime in the neutral region and other terms have their usual meanings. 


Thus, the ratio of hole current Ip and the thermionic emission-diffusion current I, given by 
eq.(l 4), is written as 


Ip 

I 


qNcNy 

A**NdT^ 




(1.17) 


Here, we have made the substitutions Lp = yDpTp and from eq.(1.15). Thus, the ratio 
{Ip /I) increases with the increase of barrier height {(j)^) and diffusion constant Dp and 
decrease of Nu, Tp and T. Yu and Snow [10] have shown that (Ip/I) is less than 0.1 for 
Schottky barriers on n-type silicon having (l)b = 0.8 V and N^, » 10^°w~^. At still large 
forward bias, signiflcant voltage drop occurs within the neutral region of the semiconductor 
and causes a drift current as well and total minority carrier current increases substantially. 
The critical current L beyond which the {Ip/ 1) increases almost linearly is given by [11] 


Ic = 


qDnAdNn 

_ 


(1.18) 


Where L is the width of quasi neutral region of the n-type semiconductor and £)„ is the 
electron diffusion constant. 


1.4 Methods of barrier height measurement 


The current through a Schottky barrier diode at a forward bias V, based on the thermionic 

emission-diffusion mechanism is given by the eq.(1.8). For the forward bias V in excess 

of 3kT/q, the plot of ln(I) vs V is linear upto a point fl-Rs) term remains ineffective 

but deviates afterwards and shows saturation. The straight line portion of ln(I) vs V plot 

extrapolated to zero bias yields ln(/s) as intercept at the ordinate. Substituting the values 

of Is and the effective Richardson constant A** in eq.(1.7), zero-bias barrier height 4>bo 

can be determined for a given diode area Ad at any temperature T. It may be mentioned 

that the value of ^bo is not very sensitive to the choice of A**, since at room temperature 

a 100% increase in A** causes a increase of only 0.018 V in (f>bo- Alternatively, one can 

make ln(Is/T^) vs 1/T plot (termed as activation energy or ArrhenQis plot) by finding 

out saturation current at various temperatures and obtain the value of (/)bo from the slope 

of the best fit straight line itself. Also, the intercept at the ordinate gives ln{A**Ad), the 

product of the electrically active area Ad and the effective Richardson constant A**. The 

Y 

principle advantage of this approach is that value of electkally active area Ad is not at all 
required for the evaluation of <f>bo- 

The barrier height can also be determined by measuring the capacitance of the Schottky 
diode as a function of reverse d.c. bias (Vr) When a small a.c voltage of a few mV is 
superimposed upon a d.c. bias, charges of opposite signs are induced on the metal surface 
and in the depletion region of the semiconductor. The relationship between the capacitance 
(C) and reverse bias Vr is given by [1] 


C = Ad^ 


CsqNo 


(1.19) 


2 (* -4’n-VR-f) 

where is the permittivity of the semiconductor, Nd is the dopant concentration, is 
the distance of the Fermi level from bottom of the conduction band and other symbols 
have their usual meaning. Thus, a plot of 1/C^ versys Vr gives a straight line with slope of 
{2/Ad(-sq^D) and an intercept on the voltage axis Vro at which 1/C^ is zero. The value 

of barrier height ^6 can then be obtained easily as 

kT 

4>b = Vro + (l>n ( 1 . 20 ) 


and (pn, tbe depth of Fermi level below the conduction band, can be computed from the 
doping concentration {Nd)- Notice that Nd can be evaluated from the slope of 1/C^ vs 
Vr plot too. 


1.5 Importance and applications 


The most important commercial application of Schottky diodes is in bipolar integrated 
circuits as clamps and to a lesser extent as load resistor substitutes, diode couplers and 
level shifters. They are also used as discrete microwave transistors, optical and nuclear 
particle detectors, etc. All the practical applications, in fact, exploit the inherent high 
speed advantages of majority carrier conduction. Metal- semiconductor rectifiers are rela- 
tively immune to the speed limitation associated with the minority carrier storage in p-n 
junction devices [12-14]. The advantages of Schottky diode over p-n junction lies in its 
fast switching speed of about 1 ns, high forward conductance per unit area, large reverse 
impedance, compact size, structure which is easily integrated and fabrication at nearly low 
temperatures The disadvantages stems from the fact that the Schottky diode is a surface 
device and consequently sensitive to surface contamination and perimeter effects 


1.6 Objective of the present work 

Extensive studies have been undertaken on different aspects of Schottky barrier diodes with 
the view to understand their characteristics better and find newer applications [3,9,15-181, 
The most important aspect of metal-semiconductor contacts is the process which de- 
termines the flow of charges over the barrier from the semiconductor to the metal and 
vice-versa. Detailed knowledge of the conduction process involved is needed to extract the 
values of the barrier parameters. The reports available in the literature provide informa- 
tion about the conduction mechanisms. But thesevary in the interpretation of the I-V 
and C-V characteristics of a variety of Schottky diodes fabricated on n- and p- type silicon 
with different metals or silicides. Chand and Kumar [19,20] have studied I-V characteris- 
tics of Pd 2 ‘S'i/silicon junctions over a wide temperature range and interpreted the results 
on the basis of thermionic emission-diffusion (TED) mechanism and assuming Gaussian 
distribution of barrier heights. Suzuki et al. [21] have studied C-V characteristics of 
Al — Si02 — SiC at different frequencies under dark and illuminated conditions. The I-V 
characteristics of Schottky diodes formed by depositing various metals such as Ir, Pt, Cr, 
Er, W, Cu, Ti and Pd on silicon have been studied in detail [22-31]. However, a few brief 
reports are available on the gold based M-S junctions [32-37]. Horvath has measured C-V 


characteristics of Au/n-GaAs Schottky contacts and interpreted them on the basis of log 
normal lateral distribution of barrier heights. Sadiq have studied Au/p-AlSb 

Schottky contacts. 

Chen et al. [34] have measured photovoltage, internal photo emission, and current- 

0 

■\jltage (I-V) and capacitance-voltage (C-V) characteristics of Au/n-Si(lll) Schottky barrier 
diodes m a temperature range of 7-300 K. Both the C-V and internal photo-emission data 
are found to yeild reliable values of barrier height and exhibit a negative tempeya.tare 
coefficient identical to that of the indirect band gap in silicon. On the other hand, photo 
voltage and I-V characteristics are shown to give lower values of barrier height due to the 
recombination current and display a strong positive temperature dependance. At room 
temperature, the barrier height obtained is 0.8V and ideality factor 77 = 1 . 1 . 

An attempt has been made here to fabricate Au/n-Si Schottky diodes and study their 
I-V and C-V characteristics. The advantage with gold is that it does not react with silicon 
and, therefore, can possibly lead to a true metal-semiconductor junction. In addition, 
simulation of I-V characteristics has been carried out on the basis of TED mechanism 
and assuming barrier inhomogeneities with a view to understand the effect of various 
barrier parameters in detail. 


Chapter 2 


Experimental Details And 
Procedures 


2.1 Fabrication of Schottky diodes 

The Au/n-Si Schottky barrier diodes were prepared on a n/n"^ silicon wafer of (111) orien- 
tation having a 8-9 // m thick epitaxial n-layer of resistivity 0 9 — cm over the underlying 

heavily doped region of resistivity 0.01 — cm. The overall thickness of the wafer was 

around 1 mm. For this, the wafers were first degreased with soap solution and then cleaned 
with organic solvents, namely trichloroethylene, acetone and methanol in succession in an 
ultrasonic cleaner and dried. Subsequently, wafer was etched in hydrofluoric acid for 1 min., 
rinsed in distilled water, and dipped in methanol until transferred to a vacuum chamber 
after drying with a hot air blower. The wafer was exposed to air for minimum possible 
time during this process. After attaining a vacuum « 10~® torr in the chamber, a thick 
film of aluminium was deposited by thermal evaporation method on the n+ side of the 
wafer and annealed in a furnace after introducing into a quartz tube at 450° C for 30 min 
under vacuum 10“® torr. Subsequently, wax was applied on aluminium deposited side 
and wafer treated with dilute HF (HF:iT20 = ITO) for removal of the silicon dioxide layer 
(usually formed on silicon), cleaned in trichloroethylene, acetone and methanol in succes- 
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sion and submerged in methanol until introduced again in the vacuum chamber for gold 
deposition on the epitaxial n-layer The quantity of gold m (gm) for thermal evaporation 
was predetermined from [38] m = 2'Kptr‘^ 

where 27r is the solid angle in which evaporation occurs at the source (i.e tungsten 
boat), p IS the density of gold (19 3 gm cmT^), t is the required thickness (i e , 400A) of 
the film and r is the distance of the wafer from the boat . Gold charge was placed in 
a tungsten boat mounted in the work chamber. The wafer was also placed suitably at a 
predetermind distance away from the boat After attaining vacuum 10~® torr as before, 
gold was deposited by passing a low voltage high amperage current through the tungsten 
boat on to the epitaxial n-layer (i.e., front side) of the silicon wafer through holes (1mm 
diameter) in the metal mask to form Schottky junctions. The wafer was later removed 
from the vacuum chamber and cut into small pieces to separate out the diodes. Each 
diode was then fixed on multipin header with silver paste keeping gold deposited portion 
upward. The circular gold portion of diodes were connected to header pins using very fine 
gold wires. Finally, one pin was connected to the body of the header by silver paste to get 
connection of the semiconductor side of the Schottky barrier diode. 


2.2 I-V and C-V measurement 

Fig 2.1 shows schematically the experimental set-up used for current- voltage (TV) mea- 
surement of the Au/n-Si^ihottky diodes. A programmable voltage source Keithley model 
230 and an autoranging picoammeter Keithley model 485 were employed for bias applica- 
tion and measurement of the current, respectively. These were in fact interfaced with a 
personal computer Zenith PC- XT 286 for ensuring automatic recording of current at volt- 
age step of O.OIV in the range 0.01-0.6 V. For achieving steady-state conditions, a delay 
of^few seconds was maintained in the subsequent measurement. Moreover, an average of 
three current readings were taken at each voltage step. The circuit diagram for the current- 
voltage measurement is shown in Fig 2.2 . The capacitance of Au/n-Si Schottky diode was 
measured using an impedance / gain-phase analyzer model HP4194A at a test frequency 
of 1 MHz under the reversed biased condition, ie., in the range of 0-4V with a step of 
O.IV. 













2.3 Phase evaluation and microstructure 


Another set of Schottky barrier diodes were fabricated in the same run or under identical 
experimental conditions but without using any mask to get large area coverage for phase 
evaluation through X-ray diffraction and microstructure examination under a scanning 
electron microscope. A Rich Seifert ISO-Debyeflex 2002 X-ray diffractometer with Cu Ka 
radiation (A = 1.5418A ) was used for identification of phase(s) and changes resulting due 
to annealing of samples at 450°C for 30 min. The microstructure was analyzed using a 
scanning electron microscope model JEOL JSM 840A in the secondary electron (SE) mode. 


2.4 Rutherford backscattering spectrometry (RBS) 

The Rutherford backscattering spectometry set-up consists of 

(a) a 2 MeV Van-de-Graaff accelerator for production and acceleration of Helium ions, 

(b) an electromagnet for dispersion of ions in terms of their mass and energy, (c) ion 
beam focussing unit, (d) separator for removal of oxygen from the AHe^ ions of the same 
energy, (e) a specimen backscattering chamber, and (f) detector assembly and analysis 
system was used. Fig 2.3 shows schematic digram of the backscattering chamber Sample 
holder plate (size 130 mm x 25 mm x 2.5 mm ) is fixed from the top in such a way 
that the ion-beam entering from one of the port is incident normally to its surface. An 
aluminium grid is placed around the sample holder and applied a negative voltage of - 
180V to suppress the emission of secondary elecrons, resulting due to ion beam irradiation, 
from the specimen surface. The ion beam current !^5-10 nA was utilized. It was 
measured and integrated using a current integrator Ortec model 401 to obtain the total 
^SLTUou'n-t of charge. The scattered ions are allowed to fall on the detector placed at 30° 
from the normal to the specimen surface. The resulting signal is fed to a charge sensitive 
preamplifier Ortec model 142, a main amplifier Ortec model 572 and biased amplifier Ortec 

\y\ SUCCESSION 

model 444^to obtain output as voltage pulses. These voltage pulses are then collected in 
the multichannel analyzer (MCA) model - ND65 operating in the pulse height analyser 
(PHA) mode. The information is finally stored in a person oj computer for display of the 
spectrum and further analysis, using suitable programs. 

The Rutherford backscattering spectrum of Au/n-Si Schottky barriers were obtained 
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Fig. 2 .3 ScLematic diagram of the Rutherford backscattering spectrometry (RBS) chamber 
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using AHe^ ions of 1.3 MeV at a scattering angle {6 ) of 30°, the detector solid angle 
(i7) of 2 57 milli-steridian and beam dose of 2;uC It is customary to plot the normalized 
backscattermg yield with respect to channel number/ ion energy. The normalized yield is 
given by Normalized Yield = (raw counts x corr) /{Q x dfl x 5E) 

where raw counts represent the actual counts recorded, corr is the correction factor to 
compensate for inaccuracies during charge integration, Q is the beam dose in /iC, dQ is the 
solid angle in milli-steridian and 5E is the calibration factor of the MCA in KeV/channel. 
The information about the depth distribution of the diffusing element in the sample and/ or 
thickness of the deposited film can be obtained from the RBS specir«-um[38]. 
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Chapter 3 

Results and discussion 


3.1 Forward I-V characteristics 

The forward current-voltage (I-V) characteristics of Au/n-Si (111) Schottky barrier 
diodes at 300 K are shown in Figs. 3.1 and 3.2. The plot is linear over several 
order of current, indicating thereby the dominance of thermionic emission- diffusion 
mechanism. At higher forward bias, ln(I) vs V plot deviates from the straight Hne 
because voltage drop across the series resistance (R^) become effective. The fitting 
procedure involves iteration and takes 7s, and Rg as adjustable parameters. A 
ln[I/(l-exp(-(iV/kT))] versus V plot is usually made to include voltage regime ‘V < 
3kT/q as well in the analysis. Needless to say that [l — exp{—qV/kT)] factor assume 
significance at low bias and cause bending of the ln(I) vs V plot (see e.g., Fig.3.1). 
The linear portion on extrapolation gives the intercept at the ordinate as In(Js). 
Also, the slope of the associated straight line yields the ideality factor. However, 
the nature of the plot usually makes the straight line portion quite subjective and 
leads to uncertainties and error in the value of ln(7s). So, a computer program 
in basic has been used to fit the experimental I-V data in thermionic emission - 
diffusion current eq.(1.8) [19]. The comi^uter program is initially run by assuming 
the ideality factor r; and series resistance Rg to be unity and zero, respectively, 
principally to obtain an approximate value of saturation cmrent (7*) and that fits 
the experimental data. Tliis Ig is then introduced into eq.(1.8) and program is run 
again to determine the values of tj and 


20 





10 







3.2 Forward n 

af ^ ^^^^^nt-Voltage (J V) ^ 

^-^^^^.^Jemperature. ^ ^ ^^^acteristics of Au/n c. , , 

~ __ ^ barrier die 




Rs in steps. Iteration continues until one finds a set of rj, and Rs values which fit the 
experimental data over as wide range of forward bias as possible with minimum deviation. 


3.1.1 Zero bias barrier height 

The zero bias barrier height is obtained by substituting the values of R, diode area 
Ad = 7.85 X 10 (corresponding to 1 mm diameter) and effective Richardson constant 
A** = 1 12 X 10^ for n-type silicon [1] m eq.(1.7). The values of (j)},o obtained 
together with the ideality factor r] and series resistance {Rs) determined by iteration for 
various diodes are listed in Table 3 1. Clearly, there is variation in (j)bo values from 0.72 to 
0 79 V. Further, the diodes fabricated can be divided into two groups : 

(i) exhibiting values of r] in the range 1.16 - 1.66 with corresponding Rg varying from 
213 - 392 n. 

(ii) having high ideality factor and series resistance in excess of 3000 9,. 

The characteristics of the latter group can possibly be attributed to the presence of an 
interfacial oxide layer either formed during the fabrication steps because of the poor vacuum 
or could not be removed completely during the etching operation. 

The high ideality factor essentially means that less current is flowing through the Schottky 
diode. The variation of the zero-bias barrier height can be explained on the basis of barrier 
inhomogeneities. As reported earlier by Werner and Guttler [40] and Chand and Kumar 
[20], the diodes invariably exhibit barrier inhomogeneities, which can be described by a 
Gaussian distribution function. Also, the value of zero-bias barrier height measured from 
I-V characteristics correspond to apparent barrier height which depend on the mean and 
standard deviation of the distribution funtion and temperature. The values of the (pbo 
listed in Table 3.1, in fact, are for diodes fabricated in different batches. These diodes 
are likely to exhibit variation in the mean and the standard deviation. The recombination 
of the electrons and holes in the depletion region is also reported to influence the total 
current of the Schottky diodes with high barrier [9,10]. Such a process normally takes 
place effectively via localized centres with energies lying near the middle of the band gap 
[8,41]. According to chen et al. [34], if the recombination mechanism contributes 30% to 
the total current at room temperature an apparent decrease of barrier height by lomer 



Tabic 3.1: Diode Parameters as obtained for Au/n-Si (111) Schottky barrier 
from I-V and C-V measurements. 


Diode Name 

IsinA) 

Rg 

(St) 

»? 

^bo 

I-V measurement 

^bo 

C-V measurement 

Diff. 
in ^bo 

sdl 

1.06 

213 

1.16 

0.79 

0.84 

0.05 

sd2 

26.51 

632 

1.26 

0.74 

— 

— 

sd3 

3.59 

397 

1.59 

0.79 

0.81 

0.02 

sd4 

4.52 

302 

1.24 

0.79 

0.84 

0.05 

sd5 

4.53 

268 

1.32 

0.79 

0.85 

0.07 

sd6 

4.52 

248 

1.66 

0.79 

0.91 

0.12 

sd7 \ 

25.01 

4949 

1.53 

0.75 

1.57 

0.82 

sd8 

25.91 

5627 

1.35 

0.74 

1.08 

0.34 

sd9 

22.03 

6392 

1.72 

0.75 

1.02 

0.27 

sdlO 

53.26 

6847 

1.86 

0.72 

1.06 

0.33 

sdll 

13.47 

3247 

1.42 

0.76 

1.05 

0.29 

sdl2 

14.21 

4313 

1.94 

0.76 

1.09 

0.33 
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should occur, if the contribution is 99 9% at 200 K, barrier height appears reduced by 
0 12V. However, recombination saturation-current (Iro) estimated using eq.(1.14) for n- 
type silicon with n, = lO^^m'S, r = 10-®s, depletion width w = O.Sfim and diode area 
■^d — 10 comes out to 2.4 x 10“^^A. This current is two orders of magnitude smaller 
than the saturation current found m the Au/n-Si (111) Schottky diodes at 300 K. It can 
therefore be said that recombination process is not all influencing the diode current and 
hence the barrier height in the present case. Instead, barrier mhomogenieties seem to 
be responsible for the observed changes in the barrier height. The real test is, of course 
possible by measuring the I-V characteristics over a wide range of temperature below 300 K 
and extracting the value of mean and standard deviation of the distribution function. Such 
studies are expected to show different values of mean barrier height and the standard 
deviation. Needless to say, the barrier height of 0.79 V found at 300 K for some schottky 
diodes is in excellent agreement with findings of chen et.al [34], whose value of (j)bo is 0.08 

V at room temperature. 

3.2 Barrier height by C-V measurement 

The barrier height of Au/n-Si(lll) Schottky diodes has also been determined by measuring 
the capacitance under the reverse bias condition. Figs. 3 3 and 3.4 show 1/C^ verses V 
plot for various diodes. Such a plot yields a straight line as per eq (1-19) with the intercept 
Vro on the voltage axis giving the barrier height as 

= VR, + <Pn + kT/q (3.1) 

where 

A kT .Nc. /o 

4>n = — ( 3 - 2 ) 

q IMd 

The barrier height evaluated for various diodes in this manner are listed in Table 3.1. 
We see that these values are invariably higher than those obtained from the I-V data. 
Moreover, there are still two groups of diodes as pointed out earlier. This result is in 
consonance with the findings of other workers [40,42]- The discrepancy has been attributed 
to the inhomogeneities in the barrier height itself. Werner and Guttler [40] solved the C- 

V equation analytically by incorporating a Gaussian distribution function for the barrier 
height inhomogeneities and shown that correspond to mean itself 
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REVERSE VOLTAGE (V) 


Fig. 3.3 1/C^ versus Voltage for Au/n-Si Schottky barrier diodes. 






3.3 X-ray diffraction 


Fig. 3 5 shows the X-ray diffraction patterns of Au/n-Si (111) Schottky barrier diode sam- 
ples prior to (a) and after (b) annealing at 450°C for 30 minutes. The diffraction peaks 
of both the patterns correspond to gold metal. Also, as deposited gold film on n-Si (111) 
exhibits a single diflPraction peak amounting to preferred orientation of grains However, on 
annealing an usual powder diffraction pattern results which means that randomly oriented 
bigger grains have emerged The 26, interplanar spacing d (A), indices of reflection (hkl) 
for both the patterns are given in Table 3.2. 


3.4 RBS and SEM studies 

In order to get insight about the nature of the interface region of Au/n-Si schottky barrier 
diodes, Rutherford backscattering spectrometry (RBS) experiments were performed in a 
manner described in section 2.4. Fig. 3 6 shows the scattering yield of ions as a function of 
energy for samples in which gold film was deposited onto n-Si substrate. There is a sharp 
peak at « 1.2MeV that corresponds to gold species. The nature of peak suggests that 
gold is concentrated in a very narrow region in terms of thickness. The RBS spectrum 
is also simulated using the code RUMP [43]. The geometry and sample structure used 
for simulation is shown in Fig 3.7. The simulated spectrum depicted by solid line in Fig. 
3.6 matches very well with the experimental one. It indicates that there is no reaction 
at the gold-silicon interface, which appears sharply defined. Fig. 3.8 shows the spectrum 
after annealing of the sample at 450°C for 30 min. Notice reduction of the peak at 1.2 
MeV and emergence of a tail at lower energies. Such an observation generally indicates 
diffusion of gold species into silicon and/or formation of gold globules during the annealing 
process. For better comparison the RBS results are presented together m Fig.3-9. Notice 
that the peak height is reduced to nearly one third. Fig. 3.10 shows typical scanning 
electron micrographs of Au/n-Si schottky diodes prior to and after annealing at 450° C for 
30 min. While Fig. 3.10a reveals no structure and possibly giving the topography of the 
silicon surface itself, Fig. 3.10b is clearly depicting formation of globules, in agreement 
with the findings of RBS studies discussed above and also reported by Campisano etal. 
[44] • Obviously agglomeration of gold occured during annealing step. Besides there is 
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Table 3.2. Crystal data of Au/-nSi Schottky barrier diodes. 
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Fig. 3. 6 Rutherford backscattering spectrum of Aii/n-Si Schottky barrier diode. 
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Fig. 3.7 Iluthc'rford backscattcriiig spectrometry (RBS) geometry and 
structm'e used for simulation. 
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O Si/Au annealed ( 450 C , 30 min) 
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Fig. 3.8 Rutherford backscattering spectrum of Au/n-Si Schottky barrier diode after an- 
nealing at 450°C' for 30 min. _ 
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Fig. 3.9 Rutherford backscattering spectra of Au/n-Si Schottky barrier diode prior to ( ) 

annealing and after annealing ( — ) at 450°C' for 30 min, i.e., Figs 3 7 & 3.8 superimposed. 
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Fig. 3 10 Scanning clccl.)'ou inicr- grapi) of Aii/n-Si Schotiky barrier diodes showing mi- 
croal.iair.l.nro : (a) i)no)- (,o and (b, afU'i annealing at 45Q°C for 30 minutes. Notice emer- 
gence of globule, s in (b). 
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possibility of some diffusion of gold taking place into the silicon substrate. Formation of 
gold globules is the reason perhaps for deterioration in the characteristics of the Schottky 
diodes. 
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Chapter 4 


Simulation studies 


4.1 Basic equations : 

The Schottky diode parameters, namely barrier height (pb and ideality factor rj are found 
to exhibit abnormal behaviour at low temperatures. Accordingly, barrier height decreases 
while ideality factor increases with decrease in temperature. The existence of barrier 
height inhomogeneities at the metal- semiconductor interface can explain satisfactorily the 
observed characteristics [44-50]. The barrier height inhomogeneities may result due to 
variation in the thickness of the deposited metal, non-uniformity of the mterfacial charges 
and/or locally defective hot regions [51]. The Gaussian distribution function that describes 
the barrier inhomogeneities adequately has the form [40,41,51,52]. 

where the term l/crV^ is called normalization constant, ^b the mean barrier height and 
a is the standard deviation. 

The total current through a Schottky diode is then given by 

-[“CXD 

I{V)= I I{<Pb,V)P{(l>b)d(l>b (4.2) 

— 00 

where I{(j)b, V) is the current for a barrier of height 4>b at a forward bias of V volts. Accord- 
ing to Chand & Kumar [20], final expression for I(V) based on the thermionic emission- 
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diffusion (TED) theory becomes 

I{V) = A,A-Thxpl^{^, - ^)]exp(|^)ll - exp{^)] (4.3) 

Further, assuming linear voltage variation of the distribution parameters, mean barrier 
height and standard deviatm a, we can write 

h = <i>bo + iV (4 4) 


and 


a = a, + (4-5) 

Where 7 = (^) and ^ = (|^) are positive voltage co-efficients. Making these substitutions 
into eq.(4.3) and neglecting the term, we get 


J(y) = hexp{^^)[l - exp{-^)] 






where 


Is = AdA**T^exp{ 


kT 


(4 6) 


(4 7) 


stands for the saturation current Here (t)ap and rjap represent the apparent barrier height 
and standard deviation and are given by 


4^ap — 4^bo 


^ap 

_L 

Pap 


2kT 


(1 - 7) + 


<yoQ^ 

kT 


(4.8) 

(4.9) 


Eq.(4.6) is of the form of current expression resulting due to TED mechanism from a 
homogeneous barrier and has (pap a.nd Pap place of <pbo V> respectively (compare, e.g., 
with eq.(1.6)). 


4.2 Simulation of diode parameters 

Simulation is mainly undertaken to study the changes caused by various parameters (ao, 
on the barrier height vs. temperature, ideality factor vs. temperature and ln{Is/T^) vs. 
1 /T characteristics. The simulated plots can be very useful in the analysis of experimen- 
tally measured I-V data. Since, in simulation, we purposely vary some parameters and 
observe their effects, ascertaining the reason for a particular feature observed experimen- 
tally becomes possible. 
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4.2. 1 Barrier height 

Eq.(4 8) IS used to simulate the apparent barrier height {(j)ap) by varying the value of the 
standard deviation (cto), i.e, increasing the extent of barrier inhomogeneities. While the 
value of (j)i)o selected is 0.80 V (known for Au/n-Si Schottky diode at 300 1^, Gq is varied 
in the range zero to 0 10 with a step of 0.02. Fig. 4.1 shows a function of T 

for various values oi a q. This clearly shows that the decrease in barrier height is caused 
by the existence of distribution (i.e, (j^ 7^ 0 ) and its effect becomes more dominant at 
low temperatures. Also, for low (Tq, barrier height initially remains constant (or decreases 
slowly ) but decreases sharply below a certain temperature. The place of inflexion moves 
towards higher temperature as the value of ctq increases. In other words, bending of the 
(pap VS T gets initiated earlier (while cooling) if inhomogeneities are severe (i.e, ao is high). 

4.2.2 Ideality factor 

Eq.(4.9) IS utilized to simulate the apparent ideality factor (ry^p) for various values of the 
voltage coefficient ^ m the range -.020 to -t-.020 while fixing Uo and 7 at 0.06 and 0.08, 
respectively. Fig. 4.2 shows variation of simulated ideality factor {'qap) with temperature 
for different values of It clearly reveals that ideality factor [rjap) increases or (decreases) 
sharply with fall in temperature if ^ is negative ( or positive). Since experiments usually 
subscribe to increase of ijap with decrease in temperature, it is obvious that ^ should es- 
sentially be negative. Also, if ^ is zero, rjap assumes a fixed value r}— 1/(1 — 7) determined 
by the voltage coefficient (7) of the barrier height at all the temperatures. If 7 is posi- 
tive 77 > 1 whereas for negative 7's, 77 < 1. Both types of cases have been observed in 
practice[20, 40,54] The latter is meaningful only for negative ^ and jy] < {\aoq^lkT\) as 
rjap for a real Schottky diode should be greater than unity [20, 40]. Further, it is obvious 
from Fig.4.2 that ideality factor begins to rise early at high temperature for greater values 
of negative ^’s. As the temperature increases, the second term of eq. (4.9) becomes less 
effective and rjap attains a saturation value of 77 = 1/(1 - 7) Fig. 4.3 illustrates the effect of 
7 on the ideality factor rjap. Clearly, the curve gets shifted towards higher rjap with increase 
in the value of 7 by a constant factor. Fig. 4 4 shows the influence of on the ideality 
factor vs temperature plots. Notice that nature of variation is similar to that found in 
Fig.4.2. This is due to that fact that ^ and ao have similar dependence on rjap in eq.(4.9) 
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Fig. 4.2 Ideality factor (i?) as a function of temperature for various valu^of voltage coeffi 
cient (^) assuming zero-bias standard deviation cio = 0.06 and voltage coefficient (y 0.08) 
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Fig 4.3 ideality factor (,) as a function of tenrpcraturc for various valuta of voltage 
coefficient M assuming eero-bias standard deviation u. = 0.06 and voltage coeffiorent 
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Fig. 4.4 Ideality factor [ t ]) as a function of temperature for various vah 
dard deviation assuming voltage cofficient 7 and £ as 0.08 and 0.01 resnf 




4.2.3 Saturation current 

Fig 4.5 shows the ln(Js/T^) vs 1000/T plot simulated using eq.(4.7) by taking ^bo — 0 80V 
and (7o between zero and 0.08. It demonstrates that for Co = 0, i.e., for a homogeneous single 
barrier, the so called Arrheneous plot is a straight line, whose slope just corresponds to 
(f)bo = 0.80V. However, as barrier inhomogeneities (i.e., standard deviation of the Gaussian 
distribution) increases the plot deviates from linearity in the progressive manner - a result 
usually observed in practice [20,40]. The activation energy is therefore not uniquely defined. 


4.2.4 Total diode current 


Most of the resistance in the Schottky barrier diode is attributed to the neutral region 
of the semiconductor (i.e., between the depletion region and back ohmic contact). The 
presence of a thin oxide layer and/or impurities are the other sources of resistance. A 
series resistance Rs is therefore associated for each diode. So, a significant voltage drop 
occurs across Rs at large forward currents. This amounts to a reduced voltage in the 
barrier region vis-a-vis actually applied to the terminals of the diode. This is accounted 
for by replacing V by (V - IRs) in eq.(4 6). Thus, 


. ,q(V-IRs),,. -q{V-IRs)^, 


(4.10) 


with Is as eq.(4.7). Since the current (I) appears on both sides of eq.(4.10), we can not 
determine its value at any forward bias by simple methods. Hence, use is made here of 
Newton-Raphson iteration formula. Accordingly, the root of equation f(x) = 0 is given by 

[55] 


X = Xo — 


fjXo) 

f'iXo) 


(4.11) 


Where f{xo) denotes the derivative of f{x) ai x = Xo- 

We can define a current function f(I) = 0 using eq.(4.10), such that 

,a(V-IRs).r. /-9(^-^^s)m _n (AA2'\ 

f{I)^I-Isexp{ ^^ ^ )]-0 i 

Introducing parameters a(= q/kT) and b{= qRJkT) into V-H-Sof eq.(4.12) and rearranging 
the terms, we get 

/(/) = I - + 4enp(c.y - 6/)(^ - 1) (4'13) 
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(4.14) 


Differentiating with respect to I, one obtains 

fif r\ 1 I ,ClV — bl . , , 1 1 

/(/)_!+ exp{—- ) - b{~ - l)I,exp{{aV - bl){— - 1)} 

Vo.p Vap ’\ap 

Let 


Fi = Isexp{ 


{aV - bl) 

Pap 


) 


and F 2 = Isexp{(aV -bl){~-l)} 

Vap 

E.^is.(4.13) and (4.14) can be written as 

f{I) = I-F^+F2 


f'{I) = 1 + — F, 

Vap 




Applying Newton-Raphson formula (4.11) for current, one writes 


(4.15) 
(4 16) 


(4.17) 

(4.18) 


I = Io 


f{Io) 


==Io 


do ~ Fi + F 2 


( 4 . 19 ) 


/W ” 1 + JJ'i - KJ - i)Ji 

where 4 is some reasonable current and /'(/o) refers to ^\i=ia Using Tailor’s expansion 
for f(I) and retaining first two terms 




(4.20) 


The sketch shown in Fig. 4.6 elaborates the logic behind the iteration method We 
take some initial guess value /<, and calculate /(/(,)• A tangent /'(/)|/=/o is drawn and 
extended to meet I-axis at I = A (say). We now calculate the value of /(A) and draw 
a tangent/' (7) l/=/i again cutting I-axis at I 2 , so that /(A) is found out. This process 
is repeated until f(I) becomes zero. Needless to say, Newton-Raphson type formula given 
by eq. (4.11) directly follows from eq. (4 20) since f(I) = 0. Fig 4.7 shows the currents 
simulated by the method enumerated above for the case of a Schottky diode having (f)f,o= 
0.8 V, ideality factor ?7ap=l 3-1' 3 temperature of 150 K for various values of series resistance. 
It reveals that ln(I) vs V plot gives a straight line with Rs= 0 but bends at high currents 
leading to saturation at any finite Rg. The initiation of bending shifts to lower current 
and linear portion shrinks with increase in the value of series resistance Rg. That means if 
current increases rapidly (or slowly) with forward bias, bending should initiate at relatively 
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CURRENT (A) 


10 


-1 


Au/n— Si Schottky contact 



FORWARD BIAS (V) 


Fig. 4.7 Simulated Current- Voltage characteristics showing effect of series resistance {Rs) 
at 150 K assuming ideality factor 77 to be unity. _ 




higher (or lower) voltage. The lower current portion still continues to be a straight line, 
however The effect of series resistance is basically to limit the current flow through the 
Schottky diode at higher forward bias. The ln(I) - V characteristics of a Schottky diode 
having a homogeneous barrier 4>bo= 0 8 V, r)ap=l and Rs = 200 0 at various temperatures 
are displayed in Fig. 4 8. The main observation are as follows: 

The plot becomes increasingly straight and shifts towards higher bias side as temperature 
is lowered. The linear portion when extrapolated to zero bias gives In(Js) i.e., saturation 
current. It decreases significantly with decrease in temperature. Further, the bending 
occurs at the same current level irrespective of temperature. This means that IR^ becomes 
effective at a particular current for a fixed Rs and therefore bending initiates progressively 
at higher bias with decrease in temperature. This information coupled with the results 
shown in Fig. 4.7 suggests that if Rs increases somehow, the bending should initiate at 
a lower current itself. Fig. 4.9 shows simulated curve of a Schottky diode having barrier 
height (fbo= 0.80 V, and Rs=0 at 200 K for various values of ideality factor. The effect 
is such that the slope of the straight line decreases with increase in the value of ideality 
factor. Also, the saturation current remains invariant as it does not depend upon riap. 
(see, e.g., eq. 4.7). The nature of plots continues to be similar at other temperatures with 
values changing appropriately, of course Further, if Rs is given some value, bending starts 
appearing in ln(I) vs V plots, as shown before in Figs. 4.7 and 4.8. Finally/ig. 4.10 gives 
a TED current of a real Schottky diode of parameters; ^bo — 0.80 V, Rs = 200 f2 at 150 K 
for various values of ideality factor as obtained by simulation using eq. 4.10 and Newton- 
Raphson iteration formula. It shows clearly that slope of linear portion of ln(I) vs V plots 
decreases with increase of ideality factor. Also, bending occurs at the same current level at 
all temperatures due to a fixed value of the series resistance {Rs)- The nature of this plot is 
very similar to I-V characteristics actually observed and discussed mChapter 3. Obviously^ 
one can proceed with the interpretation of experimental (I-V) data in the reverse way and 
determine the barrier parameters together with the level of inhomogeneities present. 


AQ 



CURRENT (A) 



FORWARD BIAS (V) 


Fig. 4.8 Simulated Current- Voltage characteristics displaying effect of temperature as- 
suming ideality factor r] and series resitance Rg to be unity and 200 fl, respectively. 


CURRENT (A) 



FORWARD BIAS (V) 


Fig. 4 9 Simulated Current- Voltage characteristics showing effect of ideality factor (rj) at 
200 K assuming series resistance Rg to be zero 
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i. r., canpub 



imi Wa a 



CURRENT (A) 



Fig. 4.10 Simulated Current- Voltage characteristics showing combined effect of ideality 

s 

factor (rj) and series fesistance at 150 K assuming series resitance to be 200^1 . 
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4.3 Algorithm of simulation 
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4.4 Computer program in C 


# include < stdio.h > 

# include< math.h > 

# define maxit 50 

double 4, lo, delta, v, a, b, fl, f2 , /* global variable */ 
float 7 ] ; 
main () 

{ 

double reduce () ; 
double I, q, k, Ad, A** ; 
float T, R/i, ^> 1)0 ; 

extern double Iq, delta, V, a, b ; 
extern float t) ; 
g = 1.60e-019 ; 
k = 1.38e-023 ; 

Ad = 7.85e - 07 ; 

A** = 1.12e + 06 ; 

4‘bo ~ 0 -^ ) 

scant 77 ) ; 

printf (” \n 77 = %/” , 77) ; 
scant k T); 

printt C\nT = %r,T) ; 
scant(”%f’, &Rs) ; 
printf C\nRs = %r,Rs) ; 
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lo = Is = Ad* A** * T'^exp{-q * ^bo/{k * T)) ; 
printf {"\nls = %le'\Is) , 
a = q/(k*T); 
b = a* Rg] 

for ( V = 0 01; V ;= 1.0; V += 0.01) { 

I = reduce 0 ; 

priiitf (”\n %4.2/ %le” , V ,1) ; 

10 = I 

} 

} 

double redxicc ( ) { 
void func (double i) ; 

int flag = 1, n = 0 ; 
double io = lo, H , c ; 
extern double delta, b, fl, £2; 
extern float r? ; 

while (flag) { 
func {io) ; 

11 = io — {io ~ /I + /2)/(l + 6 * f^/v — b* f2* {l/rj — 1)) , 
H — hn ; 

c = floor (loglO (fabs (il) ) ) ; 

delta = .001 * pow(10, c) ; 

if ( (fabs(il-i0)<= delta )11 ( n >= maxit) ) { 

flag = 0 ; 
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rctxim(il); 

} 

else 


lO = il; 

} 

} 

void fiinc(double i) 

{ 

oxtoni double Is, a, b, fl, f2; 
extern float r/ ; 
double c; 

e = (a*v - b*i); 
fl = Is * exp(c/r/); 
f2 = Is * exp(c * (1/ tta - 1) ) ; 
return ; 

} 


Note: 77 appears vn place of rjap the above program . 



Chapter 5 


Conclusion 

1. Au/ii-Si (111) Schottiky barrier diodes can be prepared successfully by depositing a gold 
layer « 600A thickness, using a thermal evaporation technique in vacuum w 10“® torr, 
onto a thoroughly cleaned n-type silicon wafer substrate. 

2. I-V characteristics of the Au/n-Si (111) Schottky diodes can be explained on the basis of 
a thermionic emission-diffusion (TED) mechanism. The ln(I)-V plots at 300 K exhibits 
linearity over a wide current range but deviates beyond 0.3 V and saturates due to the 
eflfectivcness of voltage drop occuring across the series resistance Rs of the diode. 

3. Several diode parameters viz., barrier height, ideality factor, series resistance and saturation 
current can be determined by fitting the I-V data in the TED current expression using a 
computer program. While the values of zero-bias barrier height lie between 0.72 and 0.79 
V, the ideality factor varies from 1.16 to 1.66. The high values of the ideality factor can 
be attributed to the prevailing barrier inhomogeneities at the gold-silicon interface The 
devices can be classified into two groups depending on the values of their series resistance 
(Rs) i.e., Rs in the range of 213-292 or in excess of 3000 fl * 

4. The values of barrier height obtained from the C-V characteristics are found to be invariably 
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liightti than those derived from I-V characteristics This discrepancy is possibly arising due 
to inhoinogcncities present at the metal-semiconductor junction. 

5. X-ray diffraction, SEM and RBS studres reveal that gold does not react wrth silicon even 
(hrring annealing of the Au/n-Si(lll) Schottky diodes at 450°C' for 30 min but drffuses into 
the semiconductor form globules. 

(). A cornimter program involving Neuiton-Raphson rteration method can simulate I-V char- 
acteristics of a Schottky barrier dtode using thermionic emission-diffusion (TED) current 
exprcssiorr and assuming a Gaussian distribution of barrier height for the inhomogeneities 
prrisent at the interface. The result^suggest that the abnormal decrease of zero-bias barrier 
height {(l>ap)i increase of ideality factor (rjap) and non-linearity in the activation energy plot 
of the saturation current with decrease m temperature are caused by the inhomogeneities 
present at the metal-semiconductor interface. Also, the increasing value of the ideality 
factor {rjap) decreases the slope of the straight line portion of Inl vs V plot The effect of 
series resistance is to show saturation in the I-V characteristics at higher bias. 
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